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Resumen 


Sánchez, J. A., Álvarez, T., Pacheco, J. G., Carrillo, L., €z 
González, R. A. (julio-agosto, 2016). Calidad del agua 
subterránea: acuífero sur de Quintana Roo, México. 
Tecnología y Ciencias del Agua, 7(4), 75-96. 


Este estudio se llevó a cabo con la finalidad de conocer el 
estado de la calidad del agua subterránea del acuífero 
sur de Quintana Roo, seguido de intentos por investigar 
la distribución espacial y temporal de los parámetros 
fisicoquímicos para identificar las zonas con mejor potencial 
aptas para consumo humano. Para este propósito se aplicaron 
en conjunto el índice de calidad del agua (ICA) y los sistemas 
de información geográfica (SIG). Los datos estudiados 
corresponden a las temporadas de lluvias de los años 2002 y 
2012. Los parámetros fisicoquímicos analizados fueron pH, 
temperatura (T *), sólidos totales disueltos (SDT), dureza 
total, sodio (Na*”), sulfatos (SO,?), cloruros (CI) y nitratos 
(NO,). En la interpretación de los resultados se emplearon 
estadísticos descriptivos y mapas de distribución espacio- 
temporal de concentraciones, así como el valor del ICA. 
Los resultados denotaron condiciones de heterogeneidad 
espacial en las concentraciones de los iones analizados a 
través del medio acuífero, debido a la propia geología del 
lugar; procesos, como la disolución de minerales de rocas 
carbonatadas, dolomías y evaporitas; mezcla de agua 
dulce-marina y una posible contaminación antropogénica, 
principalmente por aguas residuales domésticas y el uso de 
fertilizantes. Los constituyentes químicos que excedieron el 
límite permisible de la NOM-127-55A1-1994 fueron: SID 
(22% en 2002, y 42% en 2012); dureza total (60%); Na* (9.8%); 
CI (9.9%), y NO, (3%, 2012). El ICA demostró que la calidad 
química del agua subterránea para consumo humano es 
aceptable para la mayoría de los sitios estudiados. 


Palabras clave: aguas subterráneas, carst, Quintana Roo, 
índice de calidad del agua, límites permisibles, distribución 
espacial. 


Abstract 


Sánchez, ]. A., Álvarez, T., Pacheco, J. G., Carrillo, L., € González, 
R. A. (July-August, 2016). Groundwater Quality: Quintana Roo, 
Mexico, Southern Aquifer. Water Technology and Sciences (in 
Spanish), 7(4), 75-96. 


This study was conducted in order to know the status of the 
groundwater quality of the aquifer south of Quintana Roo, followed 
by attempts to investigate the spatial and temporal distribution of 
physicochemical parameters to identify to identify places with the 
best quality for drinking consume within the study area. This goal 
has been achieved with the combined use of the Water Quality Index 
(WQD and a geographical information system (GIS). The analyzed 
data correspond to the rainy seasons of 2002 and 2012. Physico- 
chemical parameters analyzed were pH; temperature (T *); total 
dissolved solids (TDS); total hardness; sodium (Na”); sulfate (SO,?); 
chlorides (CH), and nitrate (NO,). Descriptive statistics and maps of 
spatiotemporal distribution and ICA were used in the interpretation 
of the results. The results denoted conditions of spatial heterogeneity 
of the medium aquifer, due to processes such as mineral dissolution 
of carbonate rocks, dolomites and evaporates, the mixing of fresh- 
salt water and possible anthropogenic pollution, mainly by domestic 
sewage and fertilizer use. The chemical constituents that exceeded 
the allowable limit of NOM-127-SSA1-1994 were STD (22% in 
2002, and 42% in 2012); total hardness (60%); Na* (9.8% ); CI 
(9.9%), and NO, (3%, 2012). The ICA showed that the chemical 
quality of groundwater for human consumption is acceptable to 
most of the sites studied. 


Keywords: Groundwater, carst, Quintana Roo, water quality 
index, permissible limits, spatial distribution. 
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Introducción 


En la zona sur del estado de Quintana Roo, 
México, al igual que en toda la península de Yu- 
catán, el acuífero es la principal fuente de abas- 
tecimiento de agua para todas las actividades 
humanas; no obstante, es altamente vulnerable 
a la contaminación debido a sus características 
cársticas, que permiten infiltración y arrastre 
de contaminantes. Aunada a esta situación, las 
perspectivas de crecimiento demográfico y eco- 
nómico impulsado por la industria turística han 
traído consigo una creciente demanda de agua 
y, por consiguiente, una creciente generación 
de aguas residuales y desechos sólidos (CNA, 
1997). 

Aparentemente, las características propias de 
la península deberían proporcionar abundantes 
suministros de agua; sin embargo, la combina- 
ción de factores, como el clima y la hidrogeolo- 
gía, han generado un sistema hidrológico en el 
que el agua dulce es escasa, teniéndose inclusive 
un suministro restringido en algunas regiones 
(Ward, Weidie, €: Back, 1985). 

Dadas estas condiciones, se han empezado 
a desarrollar diversos programas de ordena- 
miento, con la firme intención de garantizar el 
abastecimiento de agua potable para zonas ur- 
banas, así como establecer las bases para que el 
agua que se extraiga para uso público en zonas 
rurales se mantenga en condiciones adecuadas 
de calidad y en cantidad suficiente como para 
satisfacer la demanda de manera oportuna y 
permanente. 

No obstante, para impulsar la gestión 
adecuada del agua es fundamental contar con 
información necesaria sobre el estado de este 
recurso, pues así será posible guiar la toma de 
decisiones y definir acciones. En este sentido, 
el conocimiento de la composición química y 
su distribución espacial a través de la calidad 
del agua es fundamental desde el punto de 
vista de su aptitud para el consumo humano, 
al igual que por la información que pueden 
aportar acerca de la dirección y extensión de 
los sistemas de flujo subterráneo (CCPY, 2012); 
más aún, puede indicar cambios estacionales y 
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temporales en la medida que sus constituyentes 
químicos estén presentes en el acuífero (Sajil € 
James, 2013). El índice de calidad del agua (ICA) 
es una de las herramientas más efectivas para 
resumir los datos de calidad de las aguas sub- 
terráneas en un formato fácilmente expresable 
y comprensible, capaz de ser utilizado por los 
responsables políticos y ciudadanos interesados 
(Lateef, 2011). Se define como una técnica de cla- 
sificación que engloba la influencia compuesta 
de cierto grupo de parámetros de calidad del 
agua en un solo valor general (Aly, Al-Omran, 
£ Alharby, 2014). 

En el mundo se han desarrollado diferentes 
ICA para aguas superficiales y subterráneas, 
entre los que destaca el de la National Sanitation 
Foundation (ICA-NSE), desarrollado por Brown, 
McClelland, Deininger y Tozer (1970) para ríos 
de Estados Unidos, y ampliamente empleado 
en diferentes estudios internacionales. Dinius 
(1987) desarrolló un ICA similar, teniendo co- 
mo valor agregado el planteamiento de rangos 
de clasificación basados en usos específicos, 
entre los que destaca el del consumo humano. 
Backman, Bodis, Lahermo, Rapant y Tarvainen 
(1998) crearon un índice para evaluar y asignar 
la cantidad de contaminación del agua subterrá- 
nea, y lo aplicaron en el suroeste de Finlandia 
y Eslovaquia central. De forma más reciente se 
han planteado ICA para la evaluación del recur- 
so hídrico a ser destinado al consumo humano, 
previo tratamiento, que incluye parámetros 
fisicoquímicos y microbiológicos relacionados 
con el nivel de riesgo sanitario presente en el 
agua. Avvannavar y Shrihari (2008) mostraron 
los resultados de su intento por desarrollar un 
índice mediante seis parámetros fisicoquímicos 
para evaluar la calidad del agua superficial del 
río Netravathi en India. Saeedi, Abessi, Sharifi 
y Meraji (2010) aplicaron ocho parámetros 
diferentes, incluyendo K*, Na', Ca”, Mg?*, 
SO,*, CI, pH y STD, como los componentes más 
importantes del agua para el desarrollo de un 
índice de calidad de aguas subterráneas en Irán. 
El trabajo de Ketata-Rokbani, Gueddari y Bou- 
hlila (2011) informa de la creación de ICA para 
las aguas subterráneas y los resultados de su 
aplicación, a fin de evaluar la calidad del agua 
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subterránea en El-Khairat, acuífero profundo de 
Túnez. 

En México, a principios de la década de 1970, 
se desarrolló un índice de calidad del agua com- 
puesto por 18 parámetros como sistema cuali- 
tativo para hacer comparaciones de los niveles 
de contaminación en diferentes áreas. Este ICA 
define el grado de contaminación existente en el 
agua a la fecha de un muestreo, expresado como 
un porcentaje de agua pura, donde 0% del ICA 
representa agua altamente contaminada y 100%, 
agua en excelentes condiciones. Este ICA sigue 
vigente en la actualidad. No obstante, para fines 
exclusivos de evaluación, se ha decidido utilizar 
tres parámetros indicadores: demanda bioquí- 
mica de oxígeno (DBO5), demanda química de 
oxígeno (DOO) y sólidos suspendidos totales 
(SST), que muestran la influencia antropogénica 
desde el punto de vista de afectación por la pre- 
sencia de centros urbanos e industriales. 

Si bien este ICA incluye gran cantidad de 
parámetros indicadores elegidos a partir de las 
estaciones de medición de la Red Nacional de 
Monitoreo de la Comisión Nacional del Agua, 
es de mencionar que los distintos pesos relati- 
vos (Wi) asignados a cada parámetro para su 
sistematización no corresponden con las carac- 
terísticas geohidrológicas de cada región; no se 
toma en cuenta la injerencia de cada parámetro 
en su entorno, como puede ser la geología del 
suelo e incluso el uso de suelo. Se sustenta más 
bien en una estandarización lineal homogénea, 
tomando como base que el recurso hídrico para 
el abastecimiento de poblaciones se obtiene de 
fuentes superficiales y no subterráneas, como es 
el caso de la península de Yucatán. 

En este estudio se propone un sistema es- 
timativo de la calidad del agua subterránea a 
través del ICA basado en la medición física de 
siete parámetros fisicoquímicos. Cada uno de 
ellos con diferentes pesos relativos (Wi), enmar- 
cados según el grado de contaminación del agua 
del acuífero y los diversos impactos sobre los 
usos de este recurso hídrico. 

Se ilustra y explica además el comportamien- 
to de siete parámetros fisicoquímicos, donde se 
pretende realizar una valoración cuantitativa 


y cualitativa de la calidad de las aguas subte- 
rráneas de manera espacial y temporal en el 
acuífero del sureste de Quintana Roo. 


Materiales y métodos 
Área de estudio 


El área de estudio forma parte de dos subregio- 
nes hidrológicas (RH): la 33-A, conocida como 
Bahía de Chetumal y Otras, y la 33-B, conocida 
como Cuencas Cerradas (figura 1). Se localiza 
geográficamente al sureste de la península de 
Yucatán, en las coordenadas 17” 49% a 20* 12 
de latitud norte, y entre los 87? 30” y 89? 27” de 
longitud oeste, cubriendo una extensión total 
de 37 896 420 km?. Colinda hacia el nororiente 
con el mar Caribe y hacia el suroriente por 
una línea que está definida por la bahía de 
Chetumal, por el límite nacional con Belice en 
un tramo de límite natural constituido por el 
río Hondo y enseguida por el límite político 
con Belice (UQROO, 2002). Esta área tiene un 
notable contraste topográfico. Hacia el norte y la 
franja costera se observa un relieve topográfico 
plano con elevaciones que van de O a 20 msnm, 
mientras que al sur-suroccidente se observa un 
relieve ondulado con elevaciones que van de 50 
a 340 msnm (Bauer-Gottwein, Gondwe, Char- 
vet, £ Marín, 2011). La temperatura promedio 
anual es de 26.15 *C, con temperaturas máximas 
entre 36 y 38 “C, y mínimas entre 12 y 14 *C 
(CCPY, 2012). De acuerdo con la clasificación 
climática de Kóppen-Geiger, la región cuenta 
con un clima aw2 (1); cálido subhúmedo con 
abundantes lluvias en verano y parte de invier- 
no (CNA, 2006). La precipitación media anual 
es de 1 365 mm, aproximadamente 600 mm más 
alta que el promedio nacional, el cual es 772 mm 
(JICA, 2004; Alcocer, Lugo, Sánchez, éz Esco- 
bar,1999). La temporada de lluvias se presenta 
en los meses de mayo a octubre, con una lámina 
media mensual de 150 a 220 mm, mientras que 
la temporada de estiaje o secas comprende los 
meses de noviembre a abril, durante los cuales 
las láminas llovidas son menores que 60 mm 
(UQROO, 2002). 
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Figura 1. Ubicación geográfica del área de estudio. La figura (a) muestra los sitios de monitoreo del año 2012; la figura (b) 


muestra los sitios de monitoreo del año 2002. El fondo en sombreado gris indica el Modelo Digital de Elevación (DEM), 
modificado de INEGI (2014). El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en 
WGS84 Zona 16 Norte (Q). 


Marco hidrogeológico 


El área de estudio se define como una unidad 
geológica cubierta casi en su totalidad de sedi- 
mentos calcáreos de origen marino del Terciario 
reciente (Butterlin € Bonet, 1960; López-Ramos, 
1972). Se sabe por Isphording (1975) que la 
geología estructural de la región de estudio está 
conformada por una secuencia de rocas de edad 
Terciaria, constituida de evaporitas intercaladas 
con rocas calizas y minerales, como dolomita, 
aragonita, yeso y halita (figura 2). Estas rocas 
carbonatadas son altamente permeables y so- 
lubles, por lo que al ser disueltas dejan escaso 
residuo, razón por la cual se presentan espesores 
de suelo muy reducidos, que a su vez resultan 
en la presencia de un sistema acuífero cárstico 
maduro con amplias fisuras, fracturas y cuevas 
(Beddows, Smart, Whitaker, € Smith, 2002; 
Hausman, 2009; JICA, 2004). 
Fisiográficamente, el área de estudio se ubica 
en la unidad “Planicie Central”, caracterizada 
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por tener una superficie plana, ligeramente 
rugosa y ondulada (López-Ramos, 1972). 

Por sus características físicas y litológicas, 
las rocas presentes hacia la carretera de Ucum- 
La Unión y González Ortega (figura 2) forman 
parte de un acuífero de calizas de características 
variadas y depósitos aluviales, cuya permeabi- 
lidad primaria es mayor en las porciones que 
contienen fracturamiento y carsticidad, y menor 
en los estratos de caliza masiva (CAPA, 2012a; 
Perry, Paytan, Pedersen, € Velázquez-Oliman, 
2009). 

La mayor parte de la región carece de un siste- 
ma de drenaje superficial y sólo en algunas áreas 
se ha formado una red hidrológica desintegrada 
no muy bien definida, que desaparece en resumi- 
deros o cuerpos de agua superficial. El acuífero se 
explota por medio de múltiples alumbramientos; 
los tipos de captación son norias, pozos someros 
y pozos profundos, que se utilizan para dife- 
rentes usos, como público urbano, agrícola y de 
servicios, entre otros (JICA, 2004). 
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Figura 2. Mapa geológico de la península de Yucatán, modificada de SGM (2007). El contorno en negro representa el área de 
estudio. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 
Norte (Q) 


La recarga del acuífero tiene lugar durante 
los meses de mayo a octubre (UQROO, 2002). 
De esta manera, gran parte de la precipita- 
ción pluvial se infiltra al subsuelo a través de 
fracturas, oquedades y conductos cársticos en 
las calizas y evaporitas (CCPY, 2012). El agua 
fluye por el subsuelo, alcanzando las costas por 
donde es descargada al mar, el resto se pierde 
por evapotranspiración. El espesor saturado 
de agua dulce decrece desde el interior hacia la 
costa (CNA, 2006, 2009). 

En el área de estudio se identifican dos re- 
giones de acuíferos: la primera corresponde a 
un acuífero compuesto por rocas carbonatadas 


del periodo del Eoceno-Paleoceno, ubicado 
hacia las porciones con mayor altitud del sur- 
suroccidente; la segunda está compuesta por 
rocas carbonatadas del periodo del Mioceno- 
Plioceno, distribuida a largo de la costa (JICA, 
2004). De acuerdo con Marín, Perry, Essaid y 
Steinich (2004), ambas regiones se pueden des- 
cribir como acuíferos de porosidad triple, donde 
el flujo del agua subterránea circula a través de 
la matriz, fracturas y un sistema de cavernas 
interconectadas con orientación NNO-SSE, casi 
perpendicular a la costa. 

En el subsuelo, el agua sigue diferentes tra- 
yectorias de flujo controladas por el desarrollo o 
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evolución del carst profundo (Velázquez, 1986). 
La dirección del flujo de agua subterránea se 
presenta de SO-NE en las zonas de mayor alti- 
tud y de O-E hacia la costa (Marín, Perry, Essaid, 
£z Steinich, 2001; Gondwe, Hong, Wdowinski, 
éz Bauer-Gottwein, 2010; Charvet, 2009; JICA, 
2004; CNA, 2001, 2007, 2009; Perry, Velazquez- 
Oliman, € Marin, 2002). La tendencia es en re- 
lación con la mayor densidad de fracturamiento 
de las formaciones geológicas que lo conforman 
(figura 3), moviéndose de las zonas de mayor 
precipitación hacia la costa, constituyendo hacia 
esta zona su área de descarga, alimentando a su 


vez al sistema Bacalar-Río Hondo, la bahía de 
Chetumal y mar Caribe (Velázquez, 1986; CNA, 
2007). 


Metodología 


Los datos fisicoquímicos de la calidad de las 
aguas subterráneas corresponden a dos perio- 
dos de muestreo (años 2002 y 2012) realizados 
durante la temporada de lluvias. La información 
del periodo 2002 forma parte del estudio de 
Aguas Subterráneas Serie Il, realizado por el 
Instituto Nacional de Estadística y Geografía 
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Figura 3. Mapa de flujo de agua subterránea a escala regional modificado de Sánchez-Sánchez, Álvarez-Legorreta, Pacheco- 


Ávila, González-Herrera y Carrillo-Bibriezca (2015). Las flechas indican las direcciones de flujo de agua subterránea según la 


referencia respectiva. El contorno en negro representa el área de estudio. El sistema de coordenadas se muestra en Universal 
Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 
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(INEGL, 2002), del cual se obtuvieron 103 datos 
fisicoquímicos cubriendo la región de estudio. 
La información del periodo 2012 fue obtenida a 
partir un muestreo realizado por CAPA (2012b), 
cubriendo un total de 120 datos fisicoquímicos. 

Según los informes de INEGI (2002) y CAPA 
(2012b), todas las muestras fueron colectadas 
de pozos de abastecimiento uniformemente 
distribuidos en el área de estudio. Todos los 
sitios contaban con su ubicación geográfica 
(coordenadas X, Y), obtenidas a partir del Sis- 
tema de Posicionamiento Global (GPS, por sus 
siglas en inglés) y representadas en el Sistema 
de Coordenadas Métricas UTM (del inglés, 
Universal Transverse Mercator). Además, estos 
informes mencionan que todas las mediciones 
de temperatura ("C), pH y conductividad eléc- 
trica (CE) se verificaron en campo, usando un 
Hydrolab Quanta G, y los análisis químicos se 
determinaron en laboratorios certificados em- 
pleando las técnicas analíticas que establecen 
las Normas Oficiales Mexicanas (cuadro 1). Las 
muestras analizadas en el laboratorio para sus 
constituyentes químicos se determinaron de la 
siguiente manera: alcalinidad total y alcalinidad 
a la fenolftaleína por titulación con ácido sulfú- 
rico 0.02N (NMX-AA-36-2001), dureza total y 
dureza de calcio por titulación con EDTA (NMX- 
AA-072-SCFI-2001), cloruros por titulación con 
nitrato de plata (NMX-AA-073-SCFI-2001), 
sólidos disueltos totales por el método gravi- 
métrico (NMX-AA-034-SCFI-2001), nitratos por 
el método de reducción con cadmio cuperizado 
(NOM-AA-079-SCFI-2001), sulfatos por el mé- 


todo turbidimétrico (NMX-AA-074-SCFI-1981), 
sodio por el método potenciométrico (ion se- 
lectivo). 


Índices de calidad del agua (ICA) 


El índice de calidad del agua fue desarrollado 
por Brown et al. (1970) y mejorado por Deininger 
(Scottish Development Department, 1975). El 
ICA se define como una técnica de clasificación 
que engloba la influencia compuesta de cierto 
grupo de parámetros de calidad del agua en un 
solo valor general (Aly et al., 2014). El método 
ICA ha sido utilizado debido a su capacidad de 
expresar de forma integral la información de la 
calidad del agua subterránea de alguna región 
determinada (Sadat-Noori, Ebrahimi, éz Liaghat, 
2014). 

El ICA fue calculado con el fin evaluar la 
calidad del agua subterránea de la región de 
estudio de acuerdo con los límites máximos 
permisibles establecidos en la Norma Oficial 
Mexicana (NOM-127-SSA1-1994). Se utilizaron 
siete parámetros químicos de calidad de agua 
(pH, SDT, dureza total, Na*, Cl, SO,? y NO, y 
el cálculo consistió de tres etapas. En la primera, 
a cada uno de los siete parámetros se les asignó 
una ponderación (1w,) acorde con su importancia 
relativa en la calidad del agua para consumo 
humano, así como de las condiciones geohidro- 
lógicas de la región de estudio (zona cárstica). 
La ponderación máxima fue de cinco y se asignó 
a sólidos totales disueltos y nitratos, en función 
de los efectos adversos en la salud humana 


Cuadro 1. Técnicas analíticas según las Normas Oficiales Mexicanas. 


Parámetro Método 


Alcalinidad total y alcalinidad a la 


Titulación con ácido sulfúrico 0.02N 


Norma oficial mexicana 


NMX-AA-36-2001 


fenolftaleína 
Dureza total y dureza de calcio Titulación con EDTA NMX-AA-072-SCEI-2001 
Cloruros Titulación con nitrato de plata NMX-AA-073-SCEI-2001 


Sólidos totales disueltos 


Método gravimétrico 


NMX-AA-034-SCEI-2001 


Nitratos 


Método de reducción con cadmio 


NOM-AA-079-SCEI-2001 


cuperizado 
Sulfatos Método turbidimétrico NMX-AA-074-SCE[-1981 
Sodio Espectrofotómetro de absorción atómica NMX-AA-051-1981 
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(Camargo € Alonso, 2007; Pacheco, Cabrera, 
€ Pérez, 2004). A los iones sodio y cloruro se 
les asignó una ponderación de cuatro, ya que 
son indicativos de mezcla con agua marina 
(intrusión salina) en zonas costeras. La ponde- 
ración de tres fue asignado a la dureza total y 
al ion sulfato; al pH se le asignó la ponderación 
mínima de uno debido a su variación mínima. 

En la segunda etapa se determinó el peso 
relativo (W;,) de cada parámetro de acuerdo 
con la ecuación del método del peso aritmético 
(Brown et al., 1970; Horton, 1965): 


ma (1) 


Donde W, es la ponderación relativa, w, es 
el peso de cada parámetro y n es el número de 
parámetros. 

En la tercera etapa se estableció una escala 
de valoración de la calidad para cada paráme- 
tro (qi). Para esto, se dividió la concentración 
resultante de cada muestra de agua analizada 
por su respectivo límite máximo permisible de 
acuerdo con NOM-127-55A1-1994 y el resultado 
fue multiplicado por 100: 


gi= =x100 (2) 


Donde q,es la calificación de la calidad, C,es 
la concentración de cada parámetro químico en 


cada muestra de agua en mg/l y S, es el valor 
estándar de agua potable en mg/l (cuadro 2). 

Finalmente, para calcular el ICA se determi- 
nó el Slí usando la ecuación (3): 


Sli=(Wi)(gi) (3) 


Donde Sli es el subíndice de cada parámetro 
y q, es la calificación de la calidad: 


ICA= Y sti (4) 


Donde ICA es el índice de calidad del agua y 
SIi es el subíndice de cada parámetro. 

El resultado final son valores clasificados en 
cinco categorías de calidad del agua: excelente 
(< 50); buena (de 50 a < 100); pobre (de 100 a 
< 200); muy pobre (200-300), y no apta para 
consumo humano (> 300). 

Se realizó un análisis estadístico descriptivo 
para comparar los valores de las medias, desvia- 
ciones estándar y varianza de cada parámetro 
(cuadro 3). Para determinar si existe diferencia 
en la calidad del agua entre las muestras de 
agua correspondientes a los años 2002 y 2012, 
se hizo la prueba de Kolmogorov-Smirnov 
debido a que los valores no se ajustaron a una 
distribución normal. Dichos análisis estadísticos 
se efectuaron a través del programa estadístico 
Statgraphics Centurion, v. XV. 


Cuadro 2. Pesos relativos de los parámetros considerados para el cálculo del índice de calidad de agua subterránea según su 


importancia para la salud humana. 


Límite máximo permisible 0 Ponderación 

Z R 2 Ponderación A 

Parámetros Unidad México (oi) relativo 
NOM-127-SSA1-1993 sl Wi 
pH Rango pH 6.5 8.5 1 0.04 
SDT mg/l 1 000 5 0.20 
Dureza total mg/l 500 3 0.12 
NO, mg/l 10 5 0.20 
Cr mg Al 250 4 0.16 
SO; mg/l 400 3 0.12 
Na* mg/l 200 4 0.16 
Total 25 1.00 
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Cuadro 3. Estadística descriptiva del agua subterránea durante la temporada de lluvias de 2002. 


Número AS Coeficiente de 
. pol e e Desviación AN 
Parámetro Unidad de Mínimo | Máximo | Media B variación Sesgo 
estándar 
muestras (%) 
pH 103 7.15 8.23 7.63 0.21 2.81 0.32 
T e 103 25.1 31.3 27.25 1.04 3.80 1.25 
SDT mg/l 103 228.2 2 877 847.15 448.51 52.94 2.2 
pa bas 103 140.61 | 2304.84 | 693.49 381.03 54.94 1.64 
CaCO, 

Na* mg/l 103 13.79 531.07 89.83 84.64 94.22 2.98 
Cr mg/l 103 197 911.07 1892 144.76 103.99 3.47 
SO, mg/l 103 4.32 1 845.79 326.31 385.52 118.15 1.82 
NO, mg/l 103 0.062 8.8 1.6827 1192 114.09 1195 


Uso de los sistemas de información 
geográfica (SIG) 


En este estudio, la modelación de los patrones 
de distribución espacial de cada parámetro 
de calidad del agua se desarrolló con base en 
técnicas geoestadísticas, usando la interpolación 
Kriging. A través de estas técnicas se pueden 
estimar los valores espaciales en aquellos sitios 
no muestreados o en el que no se ha hecho al- 
gún tipo de medición. El método Kriging se sus- 
tenta sobre modelos estadísticos que incluyen 
la autocorrelación y generación de superficies 
de predicción con elevado grado de certeza 
(Goovaerts, 1997). Kriging comprende tres 
etapas: análisis estadístico-exploratorio de los 
datos, modelo del variograma (ajuste teórico) y 
creación de la superficie de predicción (García, 
2004). Su ecuación general es la siguiente: 


Z(0) = Ya z(x) (5) 


Donde Z (0) es la predicción de la variable de 
interés, Z, en el punto X; Z(X.) es el valor obser- 
vado en el punto (base) X, coni=1,2,...N, y A, es 
un coeficiente de ponderación con ¡ =1, 2,...N. 
Estas ponderaciones se estiman considerando la 
correlación espacial y configuración de observa- 
ciones a través de los modelos del variograma 


equiparados a los variogramas experimentales 
(Castiglioni, Castellarin, £ Montanari, 2009). 

Se elaboraron los mapas de distribución 
espacial de cinco variables hidrogeoquímicas 
(dureza total, Na”, CI, SO,?, NO,;). Los mapas se 
crearon utilizando el programa ArcGIS 10.3. El 
ajuste teórico del variograma se realizó a través 
del programa GS plus, versión 9. 


Resultados y discusiones 
Análisis estadístico descriptivo 


El resumen estadístico de la composición quími- 
ca de las muestras del agua subterránea para los 
años 2002 y 2012 se tiene en los cuadros 3 y 4. 
Las muestras de agua subterránea presentaron 
valores de pH de 6.71 a 8.23, indicando que el 
agua fue moderadamente alcalina para los pe- 
riodos estudiados. No obstante, ninguna de las 
muestras analizadas estuvo fuera de los límites 
permisibles de la NOM-127-55A1-1994. 

Los STD para el periodo 2002 variaron entre 
228.2 y 2 877.0 mg/l, con una media y desvia- 
ción estándar de 847.15 + 448.51 mg/ 1. Para el 
periodo 2012, el contenido de STD varió entre 
288.6 y 7 466.85 mg/l, con una media de 1 182.04 
+ 880.26 mg/l (figura 4a). Además de exhibir 
una gran dispersión, se observó para el periodo 
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Cuadro 4. Estadística descriptiva del agua subterránea durante la temporada de lluvias de 2012 


Núm Desviación O 
Parámetro Unidad 4 Mínimo | Máximo | Media Ñ de variación | Sesgo 
muestras estándar 
(%) 
pH 120 6.71 8 7.18 0.22 3,13 0.68 
T 2 120 21.0 32 27.64 1155) 4.70 -1.09 
SDT mg/l 120 288.6 7 466.85 1182 880.26 74.47 4.1 
o 120 122.5 2160 | 658.02 | 380.72 57.86 1.19 
Caco, 
Na* mg/l 120 1.0 1 060 106.44 125.92 118.31 4.54 
Cr mg/l 120 14.0 1578 148.65 188.65 126.91 5.18 
so; mg/l 120 0.01 2 000 334.15 393.02 117.62 2.05 
NO, mg/l 120 0.002 9.36 1.47 1.98 139.74 2.14 


2002, que más de 22% de las muestras rebasó 
el límite máximo permisible establecido en la 
NOM-127-SSA1-1994 (1 000 mg/l), mientras 
que para 2012, más de 42% de las muestras 
excedió dicho límite (figura 4b). En general, las 
menores concentraciones de SDT se localizaron 
en la planicie interna del área de estudio, con 
valores entre 500 y 1 000 mg /l, incrementándose 
gradualmente hacia las costas hasta alcanzar 
concentraciones mayores a 2 800 mg/l para 2002 
y mayores a 7 400 mg/l para 2012 (figura 4b). 
La magnitud y diferencia de los valores encon- 
trados demuestran el efecto de salinización que 
da origen a aguas de mala calidad en algunas 
regiones del área de estudio, en especial en la 
zona norte. Esto coincide con lo reportado por 
CCPY (2012), en el que reportan que la gran dis- 
persión del contenido de STD entre la porción 
norte y sur del área de estudio sugiere que el 
emplazamiento de las aguas subterráneas se da 
en materiales litológicamente diferentes, unos 
más solubles que otros; sin embargo, hacia la 
porción noreste, la menor variación indica la 
presencia de materiales muy solubles. 

La temperatura (T *C) no presentó mayores 
cambios en ambos años, teniendo valores me- 
dios alrededor de 27 *C, con una desviación 
estándar de +1 aproximadamente. El orden de 
abundancia de los aniones (SO,?* > CI > NO, 
fue similar para ambos años. La amplitud 
encontrada en los valores máximos y mínimos 
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de los cationes y aniones revela heterogeneidad 
espacial del medio acuífero debido al cambio 
dinámico entre los diferentes procesos hidro- 
geoquímicos presentes en el área de estudio, 
como la disolución de minerales de rocas car- 
bonatadas, dolomías y evaporitas, y la mezcla 
de agua dulce-marina. 


Análisis geoestadístico 


La prueba Kolmogorov-Smirnof (K-S) indicó 
que los valores de dureza total, Na*, CI y SO? 
no presentaron diferencias estadísticamente 
significativas para ambos periodos, por lo que 
se agruparon para formar un solo conjunto de 
datos. Por el contrario, los valores de NO, ex- 
hibieron diferencias estadísticamente significa- 
tivas (p > 0.05). Se obtuvo la media aritmética 
para aquellos valores que coincidieron en un 
mismo sitio en ambos periodos. 

A partir de los resultados de la prueba de K-S 
se construyeron los variogramas omnidireccio- 
nales para cada variable química con ajuste a un 
modelo teórico isotrópico (cuadro 5). En ningún 
caso se consideró la existencia de anisotropía 
debido al número de datos disponibles. 

Los modelos teóricos a los que se ajustaron 
los semivariogramas experimentales fueron el 
modelo esférico, modelo exponencial y modelo 
gaussiano. De éstos, el parámetro con mejor 
coeficiente de determinación (R?) fue el ion 
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Figura 4. Concentraciones de sólidos totales disueltos (a) 2002 y (b) 2012) en el agua subterránea de la zona de estudio. Los 
puntos representan los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). 
El Datum en WGS84 Zona 16 Norte (Q) . 


Cuadro 5. Modelo teórico de variograma con mejor ajuste para los parámetros hidroquímicos. 


o Modelo Pepita Meseta Rango R? 
2002 2012 2002 | 2012 | 2002 | 2012 | 2002 2012 | 2002 | 2012 
NO, Exponencial Exponencial 0.19 0.88 1.63 4.07 | 20700 | 31200 0.67 0.69 
Dureza total Esférico 0.07 0.27 53 600 0.84 
Cl Gaussiano 0.24 0.63 JS 0.89 
SO, Exponencial 0.33 0.61 24 300 0.72 
Na Esférico 0.29 0.86 77 500 0.81 
cloruro con 89%, mientras que el ion NO, del disposición de los desechos de humanos y de 
periodo 2002 fue el que menor coeficiente de animales; infiltración de aguas residuales sin 
determinación mostró con 67%. Se utilizó el tratamiento, y en pequeña proporción la pre- 
método de interpolación Kriging ordinario, cipitación (Pacheco et al., 2004; Price, 2003). En 
para la elaboración de los mapas de distribución contraste con el ion amonio, el nitrato puede 
espacial para cada variable química. moverse con libertad a través del acuífero, por 


lo que difícilmente precipita (Hounslow, 1995). 
Mapas de distribución nitratos Diversos estudios acerca de la contaminación 
por nitratos en aguas subterráneas han mos- 


La presencia de nitratos en el agua subterránea trado que por lo general existe una relación 
puede tener diversos orígenes relacionados inversa entre las concentraciones de nitratos y 
con actividades humanas, como la agricultura; profundidad por debajo del nivel freático (Free- 
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ze €: Cherry, 1979; Ritter € Chirnside, 1984). 
Pacheco (1985) reportó para una región rural al 
norte del estado de Yucatán que la tendencia de 
las concentraciones de este contaminante en el 
agua subterránea es de naturaleza estacional, 
con las mayores concentraciones asociadas con 
los centros poblacionales. 

Las figuras 5a y 5b muestran que existió una 
diferencia en la distribución espacial de NO, 
para los periodos 2002 y 2012, como resultado 
posible de las actividades agrícolas, por lo que 
sus cantidades pueden variar de una temporada 
a otra o de un año a otro (Velasco, 2012). Todos 
los valores de NO,” del periodo 2002 y 2012 
se encontraron por debajo del límite máximo 
permisible (10 mg/l) establecido en la NOM- 
127-SSA1-1994. 

Es importante señalar que cuando el nitrato 
está en concentraciones de 10 mg/l o superiores, 
el agua es indeseable para consumo humano 
debido al posible efecto tóxico que puede tener 
sobre los niños (cianosis) (Hernández, 2008; 
Mora éz Alfaro, 1999; Pacheco é€z Cabrera, 2003). 


Dureza total 


La dureza total representa la suma de las dure- 
zas individuales debidas a los iones de calcio, 
magnesio, y es expresado como mg/l de CaCO, 
(Appelo $ Postma, 1996; Manahan, 1991). Las 
aguas duras son por lo general incrustantes, 
mientras que las blandas suelen ser agresivas 
(Lillo, 2007). Dependiendo de varios aspectos, 
como pH del agua y concentración de carbo- 
natos de calcio y magnesio, pueden producir 
procesos de corrosión o incrustación en las 
tuberías utilizadas para el abastecimiento de 
aguas domésticas (Mora éz Alfaro, 1999). 

Los valores de dureza total de las aguas 
subterráneas en el área de estudio variaron entre 
122.5 y 2 304.84 mg /l como CaCO,. Los valores 
obtenidos se agruparon en tres categorías: las 
dos primeras cumplen con lo establecido en la 
NOM-127-SSA1-1994 (< 500 mg/l) y la tercera 
presentó valores que excedieron dicho límite 
(figura 6). 
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Figura 5. Concentraciones de nitratos (a) 2002 y b) 2012) en el agua subterránea de la zona de estudio. Los puntos representan 
los sitios monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en 
WGS84 Zona 16 Norte (Q). 
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Figura 6. Concentración de dureza total en el agua subterránea de la zona de estudio. Los puntos representan los sitios 
monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 
16 Norte (Q). 


En el área de estudio, 60% de las muestras 
rebasó la norma, siendo el valor máximo en- 
contrado de 2 304.84 mg/l como CaCO,. Con 
base en la clasificación de Custodio y Llamas 
(1976), el agua subterránea del área de estudio 
se clasificó como de dura a muy dura, lo que 
pudiera explicarse por la presencia de minerales 
disueltos, como carbonatos (CaCO,), dolomita 
(CaMg (CO,),) y anhidrita (CaSO,), caracterís- 
ticos del sistema cárstico de la península de 
Yucatán (CCPY, 2012; Pacheco et al., 2004). 

Es importante señalar que el problema que 
afronta la población al utilizar aguas muy du- 
ras no sólo radica en una mala cocción de los 
alimentos, o en la generación de sarro e incrus- 
taciones en los utensilios de cocina, y en la poca 


disolución del jabón, sino además predispone a 
los pobladores a padecer ciertas enfermedades. 
Por ejemplo, en el estado de Yucatán, Medina, 
Zaidi, Real y Orozco (2002) reportaron que el 
antecedente familiar y la dureza extrema (= 
400 mg/l de CaCO,) del agua empleada para 
consumo humano son factores que predisponen 
el desarrollo de la enfermedad de urolitiasis, la 
cual se incrementa de forma notable con la edad. 


Sodio 


El sodio no se encuentra libre en la naturaleza, 
la fuente común de los niveles de sodio en el 
agua subterránea se debe, en algunos casos, a 
la presencia de rocas evaporitas, al efecto de la 
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intrusión marina, o bien a la inadecuada dis- 
posición de desechos de humanos y animales 
(Basavaraddi, Kouser, € Puttaiah, 2012; Lillo, 
2007; Pacheco et al., 2004). Los valores de las 
concentraciones del sodio de las aguas subte- 
rráneas en el área de estudio oscilaron entre 
1.0 y 1 060 mg/l. Con base en el límite máxi- 
mo permisible de 200 mg/l, establecido en la 
NOM-127-55A1-1994, se estimó que 9.8% de los 
aprovechamientos rebasó la norma (figura 7). 
Con el fin de identificar los procesos geoquí- 
micos que originan la presencia del sodio en el 
agua subterránea del área de estudio, se elaboró 
un diagrama con la relación Na*/CLr en mEq/l 


(figura 8). De manera general, se observó que 
las aguas subterráneas del área de estudio se 
encontraban ligeramente enriquecidas en sodio. 
Además, la mayoría de las muestras estaba 
por debajo de 10 y 15 mEq/l para Cl y Na', 
respectivamente. 

Algunos estudios de toxicidad del sodio en 
las aguas subterráneas y su efecto adverso en la 
salud humana revelan que la ingestión de sodio 
en altas concentraciones puede causar hiperten- 
sión, falla del corazón, enfermedades renales, 
cirrosis hepática, toxemias del embarazo, en- 
fermedad de Menier (afección del oído interno 
que causa vértigo), e hipernatremia (elevado 
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Figura 7. Concentración de sodio en el agua subterránea de la zona de estudio. Los puntos representan los sitios monitoreados. 
El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 
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Figura 8. Relación del Na+ y Cl-. Los puntos representan las muestras. 


nivel de sodio en la sangre) en los niños (Ca- 
talán, 1981). Asimismo, como cloruro de sodio 
puede causar náuseas, vómitos, inflamación del 
tracto gastrointestinal, conmutación muscular, 
convulsiones e incluso la muerte (Basavaraddi 
et al., 2012). 


Cloruros 


La mayor parte del cloruro que se encuentra en 
las aguas subterráneas suele proceder de cuatro 
fuentes naturales: 1) agua marina de sedimen- 
tación atrapada en el interior de sedimentos; 2) 
disolución de la sal común y minerales afines 
que existen en formaciones de evaporitas; 3) 
concentración por evaporación de los cloruros 
existentes en agua de lluvia, y 4) disolución de 
las partículas de material sólido existentes en la 
atmósfera (Davis € De-Wiest, 1971). De forma 
antropogénica, el aumento de cloruros puede 
ser motivado por ciertos desechos industriales, 
aguas residuales domésticas, retornos de aguas 
de riego o bien residuos de actividades en zonas 
pecuarias (Catalán, 1981; CCPY, 2012). El ion 
cloruro no forma sales de baja solubilidad, no se 
Oxida ni se reduce en aguas naturales, tampoco 
es adsorbido de modo significativo ni entra a 
formar parte de procesos bioquímicos, por lo 
que se puede considerar como una propiedad 


conservativa del agua, resultando en un buen 
trazador (Manahan, 1991). 

La distribución espacial de los cloruros 
en el agua subterránea del área de estudio 
(figura 9) mostró un rango de concentraciones 
de 14.0 a 1 578 mg/l. El 9.9% de las muestras 
excedió los 250 mg/l, límite establecido en la 
NOM-127-55A1-1994. Estos valores se ubican 
geográficamente hacia las porciones noreste y 
noroeste de la región de estudio. Esta distribu- 
ción puede estar asociada con la influencia de 
aguas marinas, aunque no se descartan la rela- 
ción con aguas derivadas de retornos de riego 
O aguas residuales derivadas de asentamientos 
humanos. 


Sulfatos 


El origen de los sulfatos está asociado principal- 
mente con la disolución de sales sulfatadas de 
yeso y anhidrita (Hounslow, 1995; Lillo, 2007). 
El sulfato está sujeto a procesos de reducción, 
en especial en presencia de bacterias y mate- 
ria orgánica; además, posee la capacidad de 
asociarse con facilidad con otros iones (Hem, 
1986). En ambientes reductores (pH < 7), la 
forma reducida estable es el ácido sulfhídrico 
(H,S), mientras que en soluciones alcalinas 
predomina el ácido sulfúrico (HS) (Lillo, 2007). 
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Figura 9. Concentración de cloruros en el agua subterránea de la zona de estudio. Los puntos representan los sitios 
monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 
16 Norte (Q). 


La concentración del ion sulfato en las aguas 
de bebida ha sido motivo de atención debido a 
su acción catártica, habiéndose encontrado un 
efecto laxante de las aguas que contienen más 
de 750 mg/l (Cabrera 1986; Catalán, 1981). En 
la península de Yucatán, los sulfatos son uno 
de los principales iones que se encuentran en 
las aguas subterráneas y además difícilmente 
alcanzan la saturación debido a que su producto 
de solubilidad es muy alto (CCPY, 2012). 

En el área de estudio, los valores del ion 
sulfato se encontraron en el rango de 0.01 a 2 000 
mg /1. Dentro de este rango, 27.8% de los apro- 
vechamientos rebasó el límite de 400 mg/l que 
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establece la NOM-127-55A1-1994. La figura 10 
muestra que la zona suroeste y la franja central 
del área de estudio son las que presentaron me- 
nores concentraciones de sulfatos; sin embargo, 
existen áreas dentro de estas zonas con concen- 
traciones superiores a la establecida en la NOM. 
No obstante, estas últimas presentan de forma 
sorprendente concentraciones bajas en cloruros; 
ello sugiere que sus características provienen en 
específico de la disolución de evaporitas (yesos / 
anhidritas) y no de la mezcla con agua de mar. 
Por otra parte, se observó un incremento gra- 
dual hacia las zonas noreste y noroeste, donde 
se llegaron a alcanzar concentraciones de hasta 
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Figura 10. Concentración de sulfatos en el agua subterránea de la zona de estudio. Los puntos representan los sitios 


monitoreados. El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 
16 Norte (Q). 


2 000 mg/] de sulfatos. Este incremento puede 
estar asociado con tres factores: 1) la presencia 
de rocas evaporitas (yesos y anhidritas), lo que 
genera un incremento sustancial en la salinidad 
del agua subterránea (Perry, Velazquez-Oliman, 
£ Wagner, 2011); 2) la descarga de las aguas de 
las zonas acuíferas de las porciones altas, y 3) la 
presencia de aguas marinas. 


Índice de calidad del agua subterránea 


Para los dos periodos analizados, se encontró 
una gran variación en los índices de calidad 
del agua. Respecto a 2002, el índice más bajo 
fue de 14 y el máximo de 174 (figura 11a). El 
valor medio del ICA es de 53 y la mediana de 


47, lo que indica que para este periodo aproxi- 
madamente la mitad de los pozos se encontró 
en la categoría de “excelente”. De los aprove- 
chamientos muestreados, 58% se ubicó en la 
categoría de “excelente”; 36%, en la categoría 
“buena”, y el 6% en la categoría “media”. En 
cuanto a 2012, el índice más bajo encontrado 
fue de 18 y el máximo de 372 (figura 11b). El 
valor medio del ICA es de 74 y la mediana de 
69. Estos valores permitieron constatar una 
degradación en la calidad del agua, pues la 
mitad de los pozos se clasificó en las categorías 
“excelente” y “buena”. De los aprovechamien- 
tos muestreados, 34% se ubicó en la categoría 
de “excelente”; 48%, en la categoría “buena”; 
15%, en la categoría “media”; 2% en la categoría 
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Figura 11. Índice de calidad del agua (a) 2002 y (b) 2012 de la zona de estudio. Los puntos representan los sitios monitoreados. 
El sistema de coordenadas se muestra en Universal Transversa de Mercator (UTM). El Datum en WGS84 Zona 16 Norte (Q). 
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“pobre”, y 1% en la categoría “extremadamente 
pobre”. Estas dos últimas categorías se ubicaron 
geográficamente en la porción central este del 
área de estudio. Tal distribución puede estar 
asociada con el crecimiento poblacional de las 
dos comunidades que se localizan en esa región, 
ya que pasó de 870 habitantes en el año 2000 a 1 
500 en 2010 (INEGI, 2014). Este crecimiento dio 
origen a la generación de más aguas residuales 
domésticas. Otra causa de dicha distribución es 
el uso de agroquímicos en actividades agrícolas 
(producción de piña) que se desarrollan en estas 
entidades. Respecto a los procesos hidrogeoquí- 
micos que pudieran complementar el elevado 
ICA encontrado, los valores concuerdan con lo 
reportado por CCPY (2012), en el que se identi- 
ficaron elevadas concentraciones de sólidos di- 
sueltos totales, sulfatos y cloruros en el pozo de 
abastecimiento de una localidad de esa región, 
con la posible influencia de aguas marinas, y la 
asociación con aguas derivadas de retornos de 
riego y aguas residuales. 
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De manera general, el mapa de distribución 
del ICA mostró un patrón espacial cualitativa- 
mente similar para ambos periodos. La calidad 
del agua subterránea disminuye desde el su- 
roeste hacia el norte, exceptuando la parte sur 
suroeste, así como algunas zonas de la parte no- 
reste del área de estudio, en donde se encontró 
la zona más segura en términos de calidad de 
agua para consumo humano. En estas regiones, 
casi todos los valores del ICA se encontraron 
en la categoría de “excelente”. La disminución 
de la calidad se asocia de manera directa con 
el aumento gradual en las concentraciones de 
sulfatos, dureza total, cloruros y sodio. 


Conclusiones 


Los constituyentes químicos que mostraron con- 
centraciones que excedieron el límite permisible 
de la NOM-127-55A1-1994 fueron sólidos tota- 
les disueltos (22% en 2002 y 42% en 2012), sodio 
(9.8%), cloruros (9.9%) y dureza total (60%). 


Sánchez et al., Calidad del agua subterránea: acuífero sur de Quintana Roo, México 


El uso de los sistemas de información geo- 
gráfica en la elaboración de los mapas de distri- 
bución espacial de las concentraciones permitió 
observar que las menores concentraciones para 
los sólidos totales disueltos, cloruros, dureza 
total y sodio se localizaron en la planicie interna 
sur del área de estudio, incrementándose sus 
concentraciones gradualmente hacia la zona 
noroeste. El mapa de las concentraciones de 
nitratos mostró una diferencia en la distribución 
espacial para los periodos 2002 y 2012. 

El ICA permitió demostrar que la calidad 
química del agua subterránea con fines de 
abastecimiento es aceptable para la mayoría 
de los sitios estudiados. No obstante, el aná- 
lisis comparativo de los periodos 2002 y 2012 
indica de manera general una degradación en la 
calidad del agua subterránea con los máximos 
valores del ICA ubicados hacia la parte noreste 
del área de estudio, que corresponde con zonas 
que presentaron un crecimiento de la población 
y de zonas agrícolas. 

Este estudio demuestra con claridad que 
el uso de métodos como el valor ICA podría 
ser de gran utilidad para evaluar la calidad de 
las aguas subterráneas en sistemas cársticos, 
proporcionando valiosas herramientas a los 
tomadores de decisiones, a fin de comprender 
de forma integral el estado de la calidad del 
agua subterránea y tener la oportunidad de 
tomar decisiones adecuadas para un mejor uso 
en el futuro. 

Por último, este estudio demuestra la necesi- 
dad de identificar zonas alternas para abastecer 
agua en aquellas poblaciones donde se reportan 
los valores más altos del ICA, pues sin duda se 
están suministrando aguas que sobrepasan los 
límites permisibles de la NOM-127-5SA1-1994. 
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